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FOREWORD 

This  report  documents  a study  conducted  at  the  Naval  Weapons  Center,  China  Lake,  Calif, 
between  January  and  July  1976.  The  work  was  conducted  under  a target  acquisition  program 
supported  by  MIPR  RA  46-75,  AMCMS  Code  675702.12.86300. 

The  Joint  Technical  Coordinating  Group  for  Munitions  Effectiveness  is  sponsoring  work  on 
surface-to-surface  target  acquisition  under  its  Joint  Munitions  Effectiveness  Manual  for  the 
Surface-to-Surface  Division.  Current  tasks  include  the  summary  of  field  test  data  from  target 
acquisition  tests,  experimentation  on  target  camouflage,  and  the  collection  of  data  on  terrain  and 
foliage  masking  (intervisibility). 

This  report  is  a handbook  for  determining  line  of  sight  in  different  types  of  terrain.  It  was 
reviewed  for  technical  accuracy  by  Ronald  A.  Erickson  of  the  Naval  Weapons  Center. 
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INTRODUCTION  AND  BACKGROUND 


This  report  is  a companion  to  NWC  TP  5908,  Line-of -Sight  Handbook. 

Its  purpose  is  to  explain  how  the  data  were  collected  and  the  computa- 
tions made  to  produce  the  data  presented  in  that  report. 

The  objective  of  this  masking  measurement  program  is  to  present 
probability  of  line  of  sight  (LOS)  as  a function  of  terrain,  range,  and 
altitude.  Preparation  of  the  handbook  included  carrying  out  a literature 
search1*2  to  determine  if  the  required  data  existed  and,  if  it  did  not, 
if  there  was  a proven  technique  that  could  be  used  to  obtain  it. 

Map  studies,  field  studies,  and  models  were  examined.  Models  were 
rejected  for  use  because  their  correlation  with  realitv  was  not  known. 
Information  obtained  from  maps  has  many  advantages,  but  also  limitations — 
mainly  that  there  is  no  good  way  to  determine  the  effects  of  vegetation 
on  LOS.  There  were  some  field  data  in  existence,  but  information  was  all 
gathered  for  particular  sites,  with  no  attempt  made  to  generalize  according 
to  types  of  terrain.  However,  the  literature  search  did  result  in  dis- 
covery of  a technique  that  could  be  used,  with  some  modification,  to 
obtain  desired  data  on  LOS  in  the  field. 

As  a result  of  the  lack  of  detailed  and  generalizable  information, 
it  was  decided  to  undertake  a measurement  program  that  would  describe  the 
LOS  characteristics  of  various  types  of  terrain.  , The  method  to  be  used 
was  an  extension  of  that  used  to  measure  masking  around  each  of  the  targets 
during  the  JTF-2  flight  trials  in  1965. 1 Those  researchers  measured  the 
elevation  angle  and  range  to  objects  surrounding  the  target  which  mask  it 
from  view  (mask  objects).  From  this  information,  the  probability  of 
having  a clear  LOS  to  the  target  from  any  range  and  altitude  was  computed. 
The  NWC  measurement  program  described  in  this  report  used  the  same  type 
of  measurements,  but  from  several  sites  in  the  same  kind  of  terrain. 
Probability  of  clear  LOS  for  each  terrain  tvpe  was  then  computed.  The 
method  is  explained  in  detail  in  this  report. 


1 Naval  Weapons  Center.  A Review  of  l, round  Target  Masking  Effects,  by  Carol  J.  Burge  and  Robert  Stohlcr. 
China  Lake,  Calif.,  NWC,  June  1974.  (NWC  TP  566K.  publication  UNCLASSIFIED.) 

" Naval  Weapons  Center.  A Review  of  Surface  to-Surface  Masking  Studies,  by  Carol  J.  Burge.  China  Lake, 
Calif.,  NWC,  June  1975.  (NWC  TP  5773,  publication  UNCLASSIFIED.) 
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THEORY  AND  DATA-GATHER1NG  PLAN 


MASK  ANGLE,  RANGE,  AND  CRITICAL  ALTITUDE 

Consider  a single  radial  extending  from  a target  or  site  at  S,  as  in 
Figure  1.  The  angle  between  the  horizontal  plane  (H)  and  a line  with 
origin  at  S which  is  high  enough  to  clear  the  tree  is  mj . That  is  the 
mask  angle  of  the  tree.  The  angle  needed  to  clear  the  first  hill  is  m2. 
An  observer  standing  at  ground  level  between  S and  the  tree  can  see  a 
target  at  S.  If  he  is  between  the  tree  and  the  first  hill  at  a range  R 
from  the  target,  S,  the  observer  must  be  at  least  as  high  as  the  value  of 
R tan  m^  in  order  to  see  S.  Between  the  first  and  second  hills,  an 
observer  must  have  an  altitude  of  at  least  R tan  m2  in  order  to  see  S. 
Similarly,  an  altitude  of  at  least  R tan  m3  is  needed  to  see  S from  a 
range  which  is  beyond  the  second  hill.  These  altitudes  necessary  to  have 
a clear  LOS  to  the  site  are  called  critical  altitudes  (CA^,  CA2 , CA3, 
and  CA4  on  Figure  2).  They  are  a function  of  the  terrain  and  of  range. 
For  any  range,  R,  critical  altitude  is  equal  to  R tan  m,  where  m is  the 
mask  angle  in  effect  at  that  range.  Whether  LOS  exists  to  the  site  from 
any  range-altitude  combination  can  be  determined  simply  by  comparing  the 
given  altitude  with  the  critical  altitude  at  the  required  range. 

Now  consider  a circle  with  circumference  at  a given  range  from  a 
site,  S,  cutting  through  the  many  radials  extending  out  from  that  site 
(Figure  2).  There  is  a critical  altitude  value,  CA,  at  each  range-radial 
intersection  (four  are  shown  in  Figure  2) . The  mean  of  these  is  the  mean 
critical  altitude  for  that  range,  with  respect  to  that  site.  Assuming 
that  the  critical  altitudes  are  normally  distributed,  about  half  the 
critical  altitudes  would  be  higher  than  the  mean.  Therefore,  if  one  were 
to  travel  the  range  circle's  circumference  at  the  mean  critical  altitude, 
one  would  expect  to  have  a clear  LOS  to  the  site  about  50%  of  the  time. 
Critical  altitudes  with  higher  probabilities  of  a clear  LOS  mav  be  found 
by  adding  standard  deviations  to  the  mean,  because  the  mean  plus  or  minus 
two  standard  deviations  contains  about  95%  of  the  values  in  a normal 
distribution.  The  mean  critical  altitude  plus  two  standard  deviations 
should  be  an  altitude  with  a probability  of  about  0.975  of  having  a clear 
LOS  to  the  site. 


PROBABILITY  OF  LOS 

Referring  again  to  the  range  circle  in  Figure  2,  the  probability  of 
having  a clear  LOS  to  S from  the  circle,  at  a given  altitude,  is  the 
ratio  of  the  number  of  radials  for  which  that  given  altitude  is  higher 
than  the  critical  altitude  to  the  total  number  of  radials.  Thus,  once 
measurements  have  been  made  of  mask  angles  and  of  range  to  mask  objects, 
probability  of  LOS  for  any  range  and  altitude  can  be  computed. 
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NUMBER  OF  SITES 

The  number  of  sites  required  to  obtain  a reasonably  good  statistical 
description  of  the  LOS  characteristics  of  a specific  kind  of  terrain  was 
unknown.  To  obtain  a guideline  for  the  number  of  sites  needed,  the 
assumption  was  made  that  probability  of  LOS  could  be  described  by  the 
binomial  model.  This  assumption  was  made  because  there  are  only  two 
possibilities,  from  any  observation  point  in  space,  with  regard  to  LOS 
to  the  site:  it  either  exists  or  it  does  not.  If  the  probability  that 

it  does  exist  is  p,  then  the  probability  that  it  does  not  is  1-p.  It 
is  further  assumed  that  each  test  of  LOS  is  independent  of  the  other 
tests.  This  obviously  is  not  true  if  the  observation  points  tested  are 
too  close  together. 

Calculations  were  made  to  get  the  approximate  number  of  LOS  measure- 
ments needed  for  a 95%  chance  that  the  probability  of  LOS  estimated  from 
the  data  would  lie  within  ±0.1  of  the  actual  probability  of  LOS.  This 
number  was  determined  to  be  between  75  and  100.  Details  of  the  computa- 
tions are  included  in  Appendix  A. 

Although  it  could  be  argued  that  each  measurement  along  a radial 
could  be  counted  toward  the  75  to  100  required,  it  was  decided  to  try 
to  work  with  between  75  and  100  separata  radials  in  each  category  of 
terrain.  To  keep  measurements  as  independent  as  possible,  only  16  radials 
per  site  would  be  measured,  at  intervals  of  approximately  22.5  deg.  Five 
or  six  sites  were  needed  for  each  terrain  type,  to  obtain  the  desired 
number  of  independent  measurements. 


TERRAIN  CLASSIFICATION 

To  determine  LOS  as  a function  of  terrain,  the  various  kinds  of 
terrain  had  to  be  separated  into  categories.  Sites  were  classified 
according  to  two  properties:  contour  and  vegetation.  Contour  ranged 

from  flat  farmland  to  sharply  rolling  hills;  vegetation  ranged  from 
scattered  low  bushes  to  dense  forests  surrounding  a small  clearing. 

Final  classification  of  a site  was  determined  on  the  basis  of  how  terrain 
looked  on  topographic  maps,  in  aerial  and  ground  photographs,  and  from 
direct  observation. 


LOCATION  OF  SITES 

Areas  of  quite  homogeneous  terrain  in  each  category  were  outlined  on 
topographic  maps  and,  within  the  areas,  sites  from  which  measurements 
would  be  made  were  tentatively  marked.  An  effort  was  made  to  locate 
each  site  in  an  area  "typical"  of  the  terrain  tvpe,  with  no  uncharacter- 
istic features.  Sites  were  selected  in  an  "average"  position  within  the 
area — that  is,  not  on  the  highest  or  lowest  ground. 
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MEASUREMENT  TECHNIQUES  AND  EQUIPMENT 

Some  time  was  spent  investigating  the  possibility  of  using  a laser 
range  finder  for  making  the  required  range  and  elevation  measurements. 

The  idea  was  reluctantly  abandoned  because  of  safety  restrictions  that 
would  have  severely  limited  the  choices  of  terrain  where  measurements 
could  be  made. 

It  was  decided  to  emplov  a standard  surveying  technique  using  two 
theodolites,  which  we  will  call  Thi  and  Th2-  The  theodolites  were  set 
up,  as  shown  in  Figure  3,  both  aimed  at  the  same  point,  P.  The  elevation 
of  the  point,  with  respect  to  Th]^,  was  measured,  using  only  Thq . 


p 


FIGURE  3.  Range  and  Mask  Angle  Measuiement. 


To  determine  the  range,  the  distance  h'  between  Th^  and  Th2  is 
measured  with  the  aid  of  a subtense  bar.  A subtense  bar  is  a bar  of 
accurately  measured  length,  X,  with  a level  and  telescope  on  it  for 
sighting  on  the  theodolite  being  used  to  make  the  measurement.  The  bar 
was  mounted  on  the  Th2  tripod,  perpendicular  to  the  line  between  Th^  and 
Th2,  as  in  Figure  4a,  and  the  angular  subtense,  a,  of  the  bar  was 
measured,  using  Thl.  B'  was  calculated  from  the  equation  below. 


2 tan  a/2 
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where 

X = length  of  subtense  bar 
a = angular  subtense  of  subtense  bar. 

The  elevation  angle,  elt,  between  theodolites  was  measured  to  correct  for 
any  difference  in  their  heights.  B = B " cos  elt,  as  shown  in  Figure  4b. 
Then 

_ X cos  elt 
2 tan  a/2  * 

With  B now  determined,  the  angles  0j  and  02  were  measured.  Then,  by  the 
Law  of  Sines, 

B sin  02 

R = sin  (0^  + 0")  ’ 

where 

R = the  range  from  Th^. 


The  theodolites  used  in  this  study  were  a Kern,  model  DMK-3,  designated 
as  the  primary  theodolite  (Th^) , and  a Wild,  Model  T-2,  designated  the 
secondary  (Th2) . The  subtense  bar  was  also  a Wild  instrument. 
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DATA-GATHERING  PROCEDURE 


Tentative  sites  from  which  to  make  LOS  measurements  were  located  on 
topographic  maps.  Precise  location  of  the  site  to  be  measured  was  done 
at  the  scene.  Table  1 gives  the  names  by  which  the  various  sites  were 
designated  for  the  study  and  their  locations,  along  with  the  number  of 
radials  along  which  measurements  were  made  for  each.  Th^  was  set  up  on 
the  spot  designated  as  the  site.  The  tripod  for  Th2  was  located  at 
least  40  m (100  ft)  away  in  a spot  with  unobscured  visibility  of  Th]  and, 
if  possible,  somewhat  uphill  from  Thp  so  that  anything  visible  to  Th ^ 
would  probably  be  visible  to  Th2. 


Thi  was  leveled  by  its  operator  while  the  subtense  bar  was  mounted 
on  the  Th2  tripod.  The  angular  subtense  of  the  bar  was  measured  by  Thj 
and  recorded.  Th2  then  was  mounted  and  leveled.  The  vertical  leveling 
bubble  in  each  theodolite  was  centered,  and  the  elevation  angle,  elt, 
between  the  two  theodolites  was  measured,  each  using  the  crosshairs 
intersection  of  the  other  scope  as  target.  If  the  measurements  of  elt 
were  not  within  20  seconds  of  each  other,  measurements  were  repeated 
until  they  were.  While  the  theodolites  were  aimed  at  each  other,  the 
azimuth  scale  on  Th^  was  set  to  0 deg  and  that  on  Th2  to  180  deg.  This 
completed  the  setup  procedure. 


To  begin  making  measurements,  Thp  was  rotated  to  the  first  azimuth 
value  (radial)  indicated  on  the  data-recording  sheet.  The  crosshairs 
were  set  on  the  skyline,  the  vertical  level  adjusted,  and  the  elevation 
read.  The  Th2  operator  looked  through  the  scope  of  Thp  at  the  position 
on  the  horizon  on  which  the  crosshairs  were  set,  went  back  to  Th2,  and 
placed  its  crosshairs  on  the  same  spot  on  the  horizon.  This  process 
sometimes  required  a few  iterations.  When  it  was  agreed  that  both  theodo- 
lites were  looking  at  precisely  the  same  spot,  the  azimuth  was  read  on  Th2. 
The  nature  of  the  mask  object  was  recorded,  i.e.,  a hill,  rock  outcrop, 
or  tree. 


For  the  next  measurement,  the  scope  of  Thj  was  lowered  to  the  next 
mask  object  down  from  the  skyline  and  the  process  repeated.  This  was  done 
for  up  to  four  mask  objects  along  the  radial,  the  closest  ones  to  the 
site  being  ignored  if  there  were  more  than  four.  The  radials  were  done 
in  four  groups  of  four  opposing  radials  at  each  site,  as  shown  in 
Figure  5.  This  was  for  two  reasons:  first,  if  all  16  radials  could  not 

be  completed  (due  to  weather  or  time  limitations) , the  sampling  would  not 
be  lopsided;  and,  second,  in  case  an  undetected,  systematic  error  developed, 
it  would  be  spread  over  all  the  data,  rather  than  deforming  one  segment. 
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TABLE  1 . List  of  Tentative  Sites  for  LOS  Measurements. 


Terrain  type 

Name  of  site 

No . of 
radials 

Geographic  location 

A.  Fairly  flat  farm- 

High  Falls 

16 

Near  Ft.  Rucker,  AL 

land;  thick  forests 

High  Bluff  1 

16 

in  distance 

Allen 

16 

Toth  1 

16 

Toth  2 

16 

B.  Fairly  smooth 

Y-l 

ii 

Near  Ridgecrest,  CA 

desert  with  little 

Y-2 

13 

vegetation 

Y-3 

15 

Y-4 

16 

Rademacher  3 

16  i 

C.  Rolling  farmland; 

David  Hendricks 

16 

Near  Ft.  Rucker,  AL 

thick  forests 

Dundee 

16 

close 

Clayhatchee  1 

16 

D.  Moderately  rough 

Wilson  Canyon  1 

12 

NWC  range,  China 

desert  and  rolling 

Wilson  Canvon  2 

10 

Lake,  CA 

hills  with  little 

Mt.  Springs  Canyon  2 

vegetation 

Rademacher  1 

19 

Near  Ridgecrest,  CA 

15 

Cameron  1 

13 

Near  Monolith,  CA 

Cameron  2 

16 

E.  Fairly  flat  farm- 

Slocomb 

16 

Near  Ft.  Rucker,  AL 

land  with  thick 

High  Bluff  2 

16 

forests  close 

Clayhatchee  2 

16 

Clayhatchee  3 

16 

Clayhatchee  4 

16 

F.  Gently  rolling 

Golden  Hills  1 

16 

Near  Tehachapi,  CA 

hills  with 

Golden'  Hills  2 

16 

scattered  trees 

Stallion  Springs  3 

16 

Stallion  Springs  4 

16 

G.  Rough  desert  with 

Wilson  Canyon  3 

16 

NWC  range,  China 

little  vegetation 

Mt.  Springs  Canyon  1 

19 

Lake,  CA 

Mt.  Springs  Canyon  3 

15 

Mt . Springs  Canyon  4 

16 

Rademacher  2 

19 

Near  Ridgecrest,  CA 

H.  Sharply  rolling 

Stallion  Springs  1 

16 

Near  Tehachapi,  CA 

hills  with  thickly 

Stallion  Springs  2 

14 

scattered  trees 

10 
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Several  problems  were  encountered  using  this  method.  One  was  the 
difficulty  of  locating,  on  the  scope  of  Th2,  a precise,  tiny  spot  seen 
as  the  crosshairs  intersection  point  in  the  scope  of  Th] . Operator 
training  in  the  process  partially  overcame  this  problem.  Leveling  the 
theodolites  was  time-consuming;  the  instruments  were  re-leveled  and  the 
0-  to  180-deg  line  reset  if  either  of  them  was  bumped  or  if  suspicious 
readings  were  obtained.  It  sometimes  happened  that  the  point  at  which 
Thi  was  aimed  was  masked  from  Th2  by  a tree  or  hill.  In  these  cases. 


’i  Th^  was  rotated  a degree  or  two  until  a mutually  visible  point  was  found. 

Of  course,  the  azimuth  reading  on  the  data  sheet  of  the  primary  theodolite 
was  changed  accordingly.  A sample  data-recording  form  and  procedure 
list  are  shown  in  Appendix  B. 


At  each  site,  12  photographs  were  taken  looking  out  from  the  site, 
starting  on  the  north  radial  and  moving  counterclockwise  in  a circle. 
Aerial  photographs  were  taken  of  the  terrain  where  the  sites  were  located. 
In  most  cases  the  individual  sites  are  shown  in  the  pictures. 
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DATA  REDUCTION  AND  RESULTS 


COMPUTER  PROGRAMS 

Two  computer  programs  for  the  UNIVAC  1108  were  written  and  used  to 
reduce  the  data,  which  consisted  of  mask  angles  and  azimuth  angles,  as 
well  as  angles  for  measuring  B (see  Figure  3) . One  program  computes  and 
plots  probability  curves  and  critical  altitudes  for  individual  sites. 

The  other  program  combines  probabilities  and  critical  altitudes  from 
individual  sites  into  summary  probability  and  critical  altitude  curves 
for  each  terrain  tvpe.  A listing  of  the  programs  and  instructions  for 
their  use  are  in  Appendix  C. 


LOS  AS  A FUNCTION  OF  TERRAIN 


The  body  of  the  Line-of-Sight  Handbook  contains  critical  altitude 
curves  and  two  sets  of  probability  curves:  one  for  altitudes  below 

1000  m and  another  for  altitudes  below  5000  m,  for  each  terrain  type. 
Aerial  and  ground  photographs  and  topographic  maps  are  included  for  each 
type  of  terrain  measured.  An  example  of  all  these  items  for  one  of  the 
terrain  types  is  shown  in  Figures  6a  and  b,  7a  and  b,  8,  9,  and  10a  and  b. 
The  handbook's  appendix  contains  the  same  information  for  each  individual 
site.  An  example  of  individual  site  data  is  shown  in  Figures  11a,  b, 
and  c,  12a,  b,  and  c,  13,  and  14. 


Figure  15  shows  the  terrain  types  included  in  this  studv,  listed  in 
order  from  that  with  the  least  masking  to  that  with  the  most.  The  ranking 
is  based  on  a comparison  of  the  probabilitv  curves  from  the  various  tvpes. 
In  an  effort  to  find  measures  that  might  correlate  with  LOS  probabilitv, 
the  average  angle  between  the  skyline  and  the  horizontal  plane  and  the 
median  range  to  the  skyline  were  computed  for  each  site.  Figure  15  shows 
the  average  of  the  average  skyline  angles  of  all  the  sites  in  each  terrain 
type.  The  standard  deviations  about  these  means  are  also  shown.  The 
average  skyline  angles  were  rank-ordered  and  the  Spearman  rank  correlation 
coefficient,  rs,  was  computed  with  the  ollowing  equation: 


E d 

n (n2  - 1) 


where 


di  = the  difference  between  ranks  for  each  terrain  tvpe, 
n = the  number  of  terrain  types. 
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For  this  study,  rs  = 0.95,  which  is  significant  at  the  0.005  level.  It 
is  not  surprising  that  average  skyline  angle  and  probabilitv  of  LOS  are 
closely  correlated,  since  the  computation  of  the  probabilitv  depended 
in  part  on  the  skyline  angle.  However,  it  does  indicate  that  the  need 
to  measure  masking  objects  below  the  skyline  should  be  reevaluated 
before  more  masking  measurements  are  made.  The  effort  may  be  worth- 
while only  for  very  low  altitudes. 

In  Figure  16,  the  average  median  range  to  the  skyline  is  shown 
for  the  eight  terrain  types.  The  standard  deviations  around  these  means 
were  huge,  so  caution  must  be  used  when  interpreting  this  plot.  It  is 
interesting  to  note  the  relatively  long  median  ranges  to  skyline  for 
desert  and  mountain  foothill  sites  compared  to  farmland  sites. 


COMPARISON  OF  FIELD  AND  MAP  DATA 

A minor  objective  of  the  measurement  program  was  to  compare  the 
field  results  with  data  obtained  from  topographic  maps.  Figure  17  (a, 
b,  and  c)  plots  probability  curves  computed  by  Erickson3  on  the  same 
graphs  with  results  from  this  study.  The  terrain  categories  were  the 
same,  but  the  actual  terrain  measured  was  different.  The  agreement 
between  the  map  and  measured  curves  is  reasonably  good,  especially  for 
the  fairly  smooth  and  the  rough  cases.  This  agreement  is  no  doubt  aided 
by  the  fact  that  vegetation  was  not  a significant  factor  for  the  desert 
terrain  types.  Probability  curves  from  map  data  have  not  been  computed 
for  other  terrain  types.  These  could  readily  be  done  if  the  need  arises. 


SUMMARY 


This  report  has  presented  detailed  information  on  how  data  were 
collected  and  calculations  made  for  TP  5908,  Line-of-Sight  Handbook . 

The  concepts  of  critical  altitude  and  of  probability  of  LOS  have  been 
explored,  along  with  classification  of  terrain  types  and  selection  of 
sites  for  making  measurements.  The  equipment  used  for  making  measurements 
and  the  techniques  used  were  described.  Data  reduction  methods  and 
examples  have  been  provided. 


''  Naval  Ordnance  Test  Station.  Empirically  Determined  Effects  of  Gross  terrain  Features  Upon  Ground 
Visibility  From  Low-Flying  Aircraft,  by  Ronald  A Erickson.  China  Lake.  Calif.  SOTS,  1?  September  1961. 
INAVWM’S  Report  7779.  NOTS  IP  2760,  publication  UNCLASSIFIED.) 


MEAN  CRIT.  ALT.,  hi  PROBABILITY  PROBABILITY 

1000  2000  3000  0 0.2  0.4  0.6  0.8  1 0 0.2  0.4  0.6 


NWC  TP  5916 


ALTITUDE, M 


RANGE,  M 

□ = 0 

O = 2000 

» = 4000 

♦ = 6000 
* = 8000 
o = 1 0500 
v = 12500 
a = 14500 


FIGURE  1 1.  An  Example  of  Data  for  an  Individual  Site,  CM  2. 
Probability  of  LOS  as  a function  of  (a)  altitude  up  to  1000  m, 
and  (b)  range;  and  (c)  mean  critical  altitude  as  a function  of 
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FIGURE  12.  An  Example  of  Data  for  an  Individual  Site,  CM  2. 
Probability  of  LOS  as  a function  of  (a)  altitude  up  to  5000  nr, 
and  (b)  range,  and  (c)  mean  critical  altitude  as  a function  of 


(a)  Aerial  view.  Sites  CM  1 and  CM 


(b)  Ground  view.  Site  CM 


FIGURE  14.  Photographs  of  Sites  CM  1 and  CM 
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FIGURE  16.  Average  Median  Range  to  Skyline. 
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(c)  Rough  desert  terrain. 


FIGURE  17.  (Contd.) 
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Appendix  A 

CALCULATION  OF  REQUIRED  NUMBER  OF  LOS  MEASUREMENTS 


The  assumption  was  made  that  probability  of  LOS  can  be  described 
by  the  binomial  model,  since  for  anv  particular  sighting  there  are  only 
two  possibilities — either  LOS  exists  or  it  does  not. 

Let  p = probability  that  LOS  exists; 

then 

1-p  = probability  that  it  does  not, 
p = estimate  of  p made  from  the  data, 

a = probability  that  p lies  within  some  tolerance,  T,  of  p, 

1 = p - T, 
u = p + T. 

Then  P[1  < p < u]  = a . 


By  the  binomial  expansion. 


u , 1 - 


n + c2 


which  can  be  simplified,  by  using  the  normal  approximation  to  the 
binomial,  to 


u,  1 «<  p ± c 


J Pd  - j), 

n 


The  normal  approximation  is  good  when  p is  close  to  0.5. 


x = c JvU-r.p)  , 


where  c is  a coefficient  obtained  from  a table  of  the  normal  distribution, 
When  solved  for  n,  the  above  equation  yields 


n = 


_ c2p(l  - p) 


t 
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It  was  decided  to  strive  for  an  estimate  of  p that  would  have  a 0.95 
probability  of  being  within  0.1  of  the  true  probability;  therefore  T 
= 0.1  and  a = 0.95.  The  following  table  was  computed  for  a p of  0.5 


n 

1 

u 

100 

0.411 

0.607 

50 

0.380 

0.657 

10 

0.274 

0.886 

It  appears  that  between  50  and  100  measurements  are  needed  to  fulfill 
the  requirements  when  p is  close  to  0.5. 

If  p = 0.9,  the  following  table  shows  that  between  36  and  50  measurements 
should  be  enough.  Since  the  number  of  measurements  required  did  not  go 
up  when  p was  raised  to  0.9,  it  was  decided  that  between  75  and  100 
measurements  for  each  terrain  type  should  be  enough. 


n 

1 

u 

50 

0.815 

0.992 

36 

0.799 

1.01 

10 

0.692 

1.14 

28 


NWC  TP  5916 


Appendix  B 

PROCEDURE  SHEET  AND  DATA  RECORDING  FORM 

MASKING  MEASUREMENT  PROCEDURE 

1.  Set  up  primary  tripod  (Kern)  at  site.  Rough  level  it.  Mount  Kern 
theodolite  on  it  and  level  it. 

2.  Set  secondary  tripod  about  67  strides  (at  least  100  ft)  from  the 
site.  Mount  the  subtense  bar  on  it  and  set  it  perpendicular  to 
line  from  the  other  tripod. 

3.  Measure  and  record  the  horizontal  angular  subtense  of  the  bar, 
alpha . 

4.  Mount  and  level  the  secondary  theodolite  on  the  secondary  tripod. 

5.  Measure  and  record  the  elevation  angle  between  theodolites  using 
the  center  of  scopes  as  targets,  elt. 

6.  Set  the  0-  to  180-deg  lines  on  the  theodolites  parallel  to  the  line 
between  tripods: 

a.  From  the  primary,  set  the  secondary  at  180  deg  (the  primary 
reads  180  deg) . 

b.  From  the  secondary,  set  the  primary  at  0 deg. 

7.  On  the  primary  theodolite,  move  the  reticle  along  the  skvline  until 
the  horizontal  circle  reading  is  80  deg.  Check  the  altitude  level, 
then  read  and  record  elevation,  El  (top  and  bottom  scale)  . 

8.  Place  the  reticle  of  the  secondary  theodolite  at  the  same  point  on 
the  skyline  as  the  primary.  Read  and  record  the  horizontal  scale, 

Azs  (middle  and  bottom) . 

9.  Lower  the  primary  theodolite  reticle  to  the  next  mask  object  down 
from  the  skyline,  keeping  the  horizontal  scale  on  80  deg.  Check 
and  read  vertical  scale. 

10.  Place  reticle  of  secondary  theodolite  on  the  same  point  and  read 
horizontal  scale. 

11.  Repeat  steps  9 and  10  for  as  manv  as  four  mask  objects  along  the 
radial.  If  there  are  more  objects,  ignore  those  closest  to  the  site. 

12.  Repeat  steps  7 through  11  for  each  azimuth  shown  on  the  data 
recording  form. 
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Appendix  C 

COMPUTER  PROGRAMS 


Both  programs  were  written  in  FORTRAN  V for  a UNIVAC  1110  computer. 


PROGRAM  PLOS 


Tha  program  "PLOS"  uses  raw  data  to  compute  ranges  to  mask  objects, 


then  critical 

altitudes  and 

probabilities 

of  LOS. 

Input  Data 

Card  1 

NT  , 

NR  , 

NHF  , 

NRP  , 

NHP 

Column 

1-5, 

6-10, 

11-15, 

16-20, 

21-25 

where 


NT  = the  number  of  sites  for  which  data  are  being  submitted, 

NR  = the  number  of  ranges  at  which  probability  is  to  be  computed, 

NHF  = the  number  of  altitudes  at  which  probability  is  to  be  computed, 

NRP:  on  the  graphs  of  probability  versus  range,  a symbol  will  be 

drawn  every  NRP  range  value, 

NHP:  on  the  graphs  of  probability  versus  altitude,  a svmbol  will 

be  drawn  every  NHP  altitude  value. 

Card  2 R^  R^  ....  R^2 

Column  1-6,  7-12,  67-72 


Card  2a  R^ 
Column  1-6, 


as  many  cards  as  needed . 


R 


NR 


R is  the  range  at  which  probability  of  LOS  is  to  be  computed. 

Six  columns  are  allowed  for  each  range  and  the  numbers  must  be  right- 
justified  . 

Card  3 HF1  HF2 HF12 

Column  1-6,  7-12,  67-72 


♦ PKECgQIfG 


r4l5S  BLANK- NOT  FILMED 
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“ “»y  cards  as  „eede(J 


JustSsd"8  aK  all««d  forta'ot  «f  LOS  is  to 


where 


Card  7 

ALP, 

Column  ?6f  ^ln  sec  oe  ELT> 

3 4~6  7-10  i3f,s  *ln  Sec  m*  NTN 

315  16~19  20-21  25_,7 

5 27  3l-36 


JO 

ar*§ular  suht- 

s.conds,  “b«„se  o,  the  subtense  ^ 

T * aI<iv»tion  angle  betu  'Srses,  minutes,  and 

seconds,  between  theodo lltes  . , 

NM  - number  of  radi  . ' d<*rees,  minutes  . 

■*  ■ ■ :Pr„:Ms 

°ard  2 AZP  rae^C  ldentification  of  fh 

D’  OMSK,  the  site. 

Column  2_|  °ef  Mln  Se  AZs 

: 14-16  £»  g:27 

Cards/radla(lXNM) ' 1,e-»  one  card  for  each 

Ch  mask  object 

A7p  J rc  measured  r 

= azln>uth  angle  fr  ’ °Ur 

r ■ --  ««. .. : :ir:r thM—  «■ 

“s  ■ «*i.uth  angle  ,r  ' and  second. 

,">J  seconds.  r°”  ae“"dary  theodolite  • ’ 

Ihera  must  be  fo„r  degrees,  WnutM 

blank  card*"?"  »'  radial.  ,f 

be  iflaer ted . hat  ~"J’  »«sk  objects  „ere 

were  not 
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Output 

The  program  prints  for  each  site: 

1.  The  distance  between  theodolites,  B. 

2.  AZP,  mask  angle,  and  range  to  mask  object,  for  each  mask  object. 

3.  Average  angle  to  the  skyline  and  average  range  to  the  skyline 
mask  object. 

4.  For  each  range: 

a.  A critical  altitude  for  every  AZP. 

b.  The  mean  critical  altitude  for  the  range. 

5.  A probability  of  LOS  table  with  probability  of  LOS  for  every 
combination  of  altitude  and  range. 

The  program  plots  for  each  site  (all  on  a single  sheet) : 

1.  Probability  of  LOS  versus  altitude — a curve  for  each  range. 

2.  Probability  of  LOS  versus  range — a curve  for  each  altitude. 

3.  Critical  altitude  versus  range — a curve  for  mean  critical 
altitude  and  a curve  for  the  mean  plus  two  standard  deviations. 

The  program  punches : 

1.  Mean  critical  altitude — one  for  each  range,  up  to  12  per  card. 

2.  Mean  critical  altitude  plus  two  standard  deviations--one  for 
each  range,  up  to  12  per  card. 

3.  A probability  table  with  a probabilitv  punched  for  each 
range — altitude  combination. 


PROGRAM  AVPL 

The  program  "AVPL"  combines  and  summarizes  the  data  from  all  the 
sites  in  a given  type  of  terrain  into  a description  of  the  terrain  type. 

Input  Data 

Card  1 NSETS 

Column  1-5 

NSETS  = the  number  of  types  of  terrain  for  which  data  is  being 
submitted  in  this  run. 
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For  each  terrain  type: 

Card  1 NT  , NR  , NHF  , NRP  , NHP  , NTN 

Column  1-5,  6-10,  11-15,  16-20,  21-25,  31-36 

These  variables  are  the  same  as  those  described  for  PLOS,  Card  1,  except 
for  the  addition  of  NTN,  NTN  is  a six-character  alphanumeric  identifier 
of  the  terrain  type. 

Card  2 R2  

Column  1-6  7-12 

Range  cards — same  as  for  PLOS. 

Card  3 HF1  HF^ 

Altitude  cards — same  as  for  PLOS. 

For  each  site  (punched  by  PLOS): 

Card  1 Mean  critical  altitude  . . . 

NR  of  them,  up  to  12  per  card. 

Card  2 Mean  critical  altitude  plus  two  standard  deviations 
NR  of  them,  up  to  12  per  card. 

PattI  ^ P P P 

LOS  11  LOS  12  LOS  1 NHF 

PLOS  21  PL0S  NR, NHF 

PL0S  = probability  of  LOS. 

Output 

AVPL  prints  the  average  probability  table  for  each  terrain  type. 
AVPL  plots: 

1.  The  average  probability  versus  altitude  with  one  curve  for 
each  range. 

2.  The  average  probability  versus  range  with  one  curve  for  each 
altitude . 

3.  Critical  altitude  versus  range — one  curve  of  the  mean  critical 
altitude  and  another  for  the  mean  plus  two  standard  deviations. 
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1 » 

2 * 
J* 
4 * 
5* 

b* 

7* 

3* 

9* 

13* 
11* 
12* 
13* 
’4* 
15* 
16* 
17* 
13* 
19* 
23* 
21  * 
22* 
23  * 

2 4 * 

25* 
26* 
27. 
2 a * 
29* 
33* 
11* 
32* 
33* 
34* 
35* 
36* 
37* 


Listing  of  the  Computer  Program  PLOS 


c 

c 

c 


c c mput es  plos  as  function  of  range  ano  altitude,  input  data  are  mask 

ANGLES  AROUND  A TARGET. 

NMrNUMBER  OF  MASK  ANGLES  PER  TARGET 

NR  ^NUMBER  OF  RANGES,  NF:NUMbER  OF  FLIGHT  ALTITUDES,  NTPNU  MB  ER  OF  TARGETS 
\P  3N UMBEP  OF  PROBABILITIES  XMRXNH  F 

DIMENSION  EMSK(2G, 4), R0(2C, 41, Rt 32  I, HC(J2,2D»,HF{22», PR 06 (32, 22) 
DIMENSION  P(32I,THETA(2u),LA0ELP(6) 

DIMENSION  ALP(3).ELT(3),A2P(72,4),CMSK(4,3),AZS(4,3),A2SC(72,4) 

DIMENSION  AVEC  132 ) , VAR  t 32) ,CA2  SDt  32) 

DIMENSION  IPAK1I210) ,1 PAK 2t 213 ) 

DIMENSION  RLA6  ( 2 ) ,hLAP  (2  I 
FORMAT  (515) 

FORMAT  (12FS.2) 

FORMAT  ( 1HJ  *PLCS  PROGRAM  SITE  ',A6> 

FORMAT  ( 1H3  7hR  A NGE  rF 7 .3, 6 H“E T ER S ) 

FORMAT  ( IHJ1SHCRITICAL  ALTITUDES) 

FORMAT  (1HJ12FIC.2) 

FORMAt  (1H323HMEAN  CRITICAL  A L TI TUOE :F 7 . 3 , 6HMET ER S ) 


1DC 
1D2 
233 
2C1 
202 
2D  3 
2 34 

235 

236 

237 
233 
330 
315 
32C 
325 


FORMAT  ( 1H 1 *P°0BA5ILI TV  CF  LC S TA3LE  SITl  • , AS  ) 
FOR“AT  ( 1HU9HALT ITUOEP 13F 13 .3) 

FORMAT! 1 hO*HR  A NGE ) 

FORMAT( lh3F6,C,4  XlDF 10,3) 

FORMAT ( F 1 .3  ) 

FORMAT  ( ’PROBABILITY  OF  LOS,  SITE  •,46.'*') 
FORMAT! ’RT’.Fb.D.’S’ ) 

FORMAT! ’ALT:’,F6.3. ’S’ > 

theta ( 1 ) -3 


SB  ARr  2, 

DC  3K I 2 , 7 3 
k 1 : K - 1 

3 THETA(K)rTHETS!Kl)»l3.C*.2C982 

CALL  FR8  3IC ( • 3175  C.  BURGE  PH.  SlbT’l 
READ  I 5.1C3 )NT , NR , NHF, NRP, NHP 
REA0(5,132)<R(K),K-t,NR) 
REA0(5,132)(HF(K)  ,K  = 1,  f.HF  ) 

N“1  :N  M ♦ 1 
D09KA  1 ,NR 
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36  * 

AVEC(K>=0 

39* 

9 

V AR  (K  )rO 

4Q* 

KS=0 

41* 

004  k A 1 , NR fNRP 

42* 

KS=KS  *1 

43  * 

ENCODE! 13, 320, 6LA3) RIK ) 

44* 

CALL  HEIGHTI3.08) 

45  * 

4 

CALL  LINES (RLAb, IPAK l.KSI 

46* 

KSAO 

47* 

OOSKr 1 ,NHF ,NHP 

48  * 

K S AKS  *1 

49* 

ENC00E(11,32S.NLABIHF(K) 

5C  * 

5 

CALL  LINESIHLA6.IPAK2.KS) 

51* 

D095LT=1,NT 

52* 

peach 5, 104)  (ALPI I)  ,1=1 ,3)  , IELT  II)  ,1:1, 3)  .NH.NTN 

53* 

104 

FORmaTIIF3.0,3X,3F3.C,Ix,I3,3X,Ao) 

54  * 

ALPH=(ALP<i)*ALPt2)/6C»ALPl3>/ : 600 ) * . 0 1 745 32 9 / 2 

55  * 

ELTHriELT(i)*ELT(2)/b3*ELT(3>/  3 b JO  >*.01745329 

5o* 

3=S3AR/  <2*  ( SIN!  4LPH  ) /COSIA  LPH  ) )»COS(  ELTH  ) ) 

57* 

WRITE lb .200 ) NTS 

58* 

WRITE  lb  .499  ) 3 

59* 

499 

FORMAT! IhO‘8  IS  ’ ,F  7 . 2 , ’ME  TER  s • ) 

6C* 

WRITE  16 . 496 ) 

6 1 * 

496 

FORMAT! 1H0’*ZIMUTH  MASK  ANCLE  RANGE’) 

62* 

WRITE  16 .497  ) 

63* 

497 

FORMA  7 1 1 h • (DEG)  IOEG)  (METERS) 

64  * 

AVEM=3 

65* 

N M R -N  M 

66  * 

AVE  R = 0 

67* 

0010K  rl ,NM 

58  * 

0081-4, 1 

69* 

REAOl E , 103 ) AZP ( K ,1) , (0  MSK  ! I ,J  ) ,J  = I ,3 ) , 1 A2S (I , J ) 

70  * 

103 

FORMAT  (2X  ,F3.C,2X,JF3.u,2X  , 3F  3 ,C  , 1 OX  ) 

71* 

AZS(I,3)=2*AZS( 1,7) 

72* 

TST=2*AZS!I,3) 

73  * 

IFCTST.LT.60IG0  TO  fc 

74* 

AZS (I ,Z> AAZS ! I,2)*I . 

75* 

A 2S  (1  ,3  ) =AZS!  1,3  )-b0. 

76* 

6 

£«SK (K, I )Z( (OMSK (I, ’ )-93>* OMSK  (I ,2 )/dO«OMSK 11 ,3 

77* 

IF (EMSK (K, I ) ,LT .G)EMSK (K , I )T0 

78* 

AZSC(K,I)-»AZSII,l)*AZSII,Zt/6C‘AZS(I,3)/360J)* 

79* 

AZPIK,I)-AZP(K,I)*.2I7m53Z9 

60  * 

ANGlr AZSC ( K , 1 ) 

51* 

ANG2AAZP (K , II-AZSCl K,I ) 

52* 

R0(K,I)=ABS(B*SIMANGI)/SIN(ANG2)) 

53  * 

AZP(K,I):AZP(k,I)/.017453  29 

84* 

WRITE  16 , 5C  J ) AZP  (K,l  ) ,EMSK  IK  ,1  ) ,RO(K,  I) 

55* 

50  J 

FORMAT (2X ,F4.3,8X,F5,2,dX , F6  .3  ) 

06* 

E"SK(K,I)rE“SK(K, II*. 017453 29 

87* 

6 

CONTINUE 

c8  * 

A VE“A  AVEM*EMSK  ( K ,l|  ) 

IAGNOSTIC*  THE  TEST  FOR  ECUALITY  BETWEEN  NON -I  N TEG  EP  S MAY  NOT 
39*  IFIRD IK ,4) .£C.U>N"RrNM-l 

90*  AVE°= AV ERwRO t K ,4 ) 


• ) 

,0-1.3) 


>/ 363u I 
.0174532° 
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914 

1G 

CONTINUE 

92* 

AV EMT AV EM/ (NM*. 21745329) 

93* 

AVER=AVER/NM 

94  * 

WRITE (6 ,531 ) AVEM, AVER 

95  * 

331 

F0R“AT(1H0’AVE  SKYLINE  ANGLEa’ 

96  * 

CALL  OPNPLT 

97* 

□0  22K; 1 ,NR 

96* 

AVErO 

99* 

SSQ-3 

123* 

OClsIrl ,NM 

1-1* 

TSTrG 

1 22* 

D015JA1 ,4 

IACNCSTIC* 

ThE  TEST  FOR  ECjALITY  BETWEEN  NON 

12  3* 

IFITST.EO.l.JGO  TO  15 

IAGNCSTIC* 

ThE  TEST  FOR  EQUALITY  BETWEEN  NON 

124* 

IF (E“SK ( I ,j) ,EQ .j)GO  TO  15 

125* 

IF(P ( K 1 .LE. RD ( I , o) ) GO  TO  11 

1-6* 

IF(J.E0.4)30  TO  13 

127* 

GO  TO  15 

1‘28* 

11 

I F ( J . GT .1)30  TO  12 

109* 

ANGAO 

1 1G  * 

GO  TO  14 

111* 

12 

JJ  = J-  1 

112* 

ANG=E“SK ( I , JJ> 

113* 

GO  TO  14 

114* 

13 

ANG=E"SK(I,JI 

115* 

14 

HC (K , I) =R ( K ) *T  AN( ANG  > 

116* 

TSTrl .3 

117* 

15 

CONTINUE 

118* 

AVE9AVE  *HC ( K , 1 1 

119* 

16 

SSa9SSQ*HCtK,I)»*2. 

1 2C  » 

AVEC(K) TAVE/NM 

121* 

VAR(K)r(SSL-(AvE«»2.)/NM)/(NM- 

122* 

CA2SD(K):avEC(k)*2*(VAR|k) )**C 

123* 

002CI , N H F 

124* 

PPOB ( K , I ) :Q 

125* 

0018LA1 , NM 

126* 

18 

IF(HF(I1.GE.HC(K,L))PRCE(K,I). 

127* 

20 

PR09(K,I)rPPCB(K,I)/NM 

128* 

WRITE (6  ,220 ) NT  h 

129* 

WRITE  (6,201)R(K) 

1 32* 

WRITE  (6,202) 

1 21* 

L2rC 

132* 

LIN=INT (NM/ 12 ) 

1 23* 

0C25LA1 , LIN 

134* 

L 1 9L2  *1 

135* 

L2  9L« 12 

1 3o* 

25 

WRITE  (6  ,2C3)  (HC  (K,I  ) ,I=L1,L2) 

137* 

IF (L2.E0.NM)G0  TO  30 

126* 

LI  =L2*1 

1 29* 

WRITE (6 ,202) (HClK.I ) ,ITL1,NM) 

142* 

20 

WRITE (6,2C4)AVEC(K) 

141* 

NA:l 

142* 

CALL  PHYSOR ( 3.2, 7 ,C ) 

• ,F6. 1 . ’meters’ ) 
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CALL  SCLPICIZ.5) 

ENC0DEI35,315,LA6ELP)NTN 
CALL  INTAX.. 

CALL  TITLE!  • * . 1 , ’ A LT I TUOE  , Mi • , ICO , • PR  OB  A? IL ITr l • , 1 CJ . 2 . 5 , 1 . 75 ) 
CALL  MSSSAGI LAEELP, 100.-J.25, 2 .75) 

CALL  CRAF(0.,5UC.,5CCQ.,J.,C.2,1.0) 

005  3K :l  ,NR  , NRP 
0052J -1 , NH  F 

52  P<J)=PR0E1K,J) 

53  CALL  CUR VE ( HF , P, NHF , 2) 

NRLrINT (NR/NRP ) 

CALL  L£GEN0<IPAKI,NRL,2.73,-.2£> 

CALL  ENDGR(O) 

.RITE  (S,205)NTN 
NHF 1- INT INHF/lu ) 

L2=0 

IF!NHF.LT.1C)G0  TO  56 
CC55 I -1 , NHF 1 
L1=L2*1 
L2=’0*I 

WRITE  16  , 226) IHFIL)*L=L1,L2) 

WRITE  (6  ,207  ) 

0055K Al ,NP 

55  .PITEC6.2C5)  R ( « ) , ( PRO  b 1 K , j ) , J 1 , L2  ) 

56  NH F2- NH F - 1 0 * NHF 1 

IF CNHF2 .EL .0)  GO  TO  59 
L 1 RL2  *1 
L 2 ~L2  ♦NHF2 

•RITE <6, 205)1 HF(L),L=L1,L2> 

WRITE  16,2^7) 

DC59K =1 , NR 

56  WRITE (6 ,203 )R(K) , (PRCB (K,J) ,J=L1,L2) 

PUNCH  602 , ( A WFC I K ) , K: 1 ,NR ) 

PUNCH  6C0, (CA250(K) ,Xri,NR) 

600  FCRHAT ( 12F6.0) 

PUNCH  601,  ( (PROS ( K , J ) ,jrl,NHF ) ,K  = 1 ,NR  ) 

601  F0RMAT(12F6.3) 

59  CALL  PHYSOR f 3. ,4 .C) 

CALL  TITLE!’  ’ , 1 , ’ R A NG E , Mi  ’ , 1 r C , ’ PRO B A 6 I Li T Y i • , 100 . 2 . 5 , l . 75 ) 

CALL  GRAF!  3.  ,22  00  .,  loOGO.  , U«,  0 .2,1 .3  I 
CALL  INTAXS 
C 06  5 J - 1 .NHF ,NHP 
C060K Zl , NR 

60  p IK  )rPROb  ( K , J) 

65  CALL  CURVE(R,P»NR,4> 

NHL-INT (NHF/NHP I 

call  LEGENU1  IP  A K 2, NHL, 2.75,-C. 15) 

CALL  ENOGRIOI 
CALL  PHYSOR 1 3. . 1 .5  I 

CALL  TITLE!’  ’ , 1 , ’RANGE , M S ’ , 10  C, ’CRI TI CAL  AL T I TU DE , Mi • , 1 00 . 2 . 5 . 1 . 7 
X5  ) 

CALL  5RAFIC.,  7C0J.,  160 GC . , 0 . , 2 CO J. , 6 CO  0. ) 

CALL  INTAXS 

CALL  CURVE 1R.AWEC, NR, NRP) 

CALL  CURVEIR,CA2S0,NR,NRP) 

CALL  ENOoR I C I 
CALL  ENDPL(O) 

95  CONTINUE 

CALL  00N2PL 
END 


0 OF  COMPILATION: 
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Listing  of  the  Computer  Program  AVPL 


C PLOTS  A VE  PROG  LOS  F R OH  PROB  TABLES 

DIMENSION  AVEC(22),AVEwtJ2),CA2SD(32),CA2(32) 

D INENSION  R(32),MF(22),PR08(32,22),PROBA(32,22>,P(32),LABELP(4> 
DIMENSION  I P Ak  1 (210) , I P A K 2 (210)  , R LAB (2) f hLA 3 ( 2) 

READt  5,  10ONSETS 
ICO  F0RHATU5) 

CALL  FRSOIDt'  3175  C .BURSE  PH.  316T') 

D 09  5L  T* 1 . NS  ET  S 

READC5,101)NT,NR,NHF,NRP,NHP,NTN 
1C1  F0RMAT(5I5,5X,A6) 

READt  5, 102) (R (K) , K=1 ,NR) 

1C2  F ORMAT ( 1 2 F6 .0 ) 

READ(5,102)(HF(K>,K=1,NHF) 

C PACKS  LEB  ENDS 

< s=c 

D04K=1,NR  ,NRP 
K S = K S + 1 

ENCODE (10 ,320  ,RLAB) R (K) 

4 CALL  LINE  St  RL»B  , I PAK  1 ,-„S) 

320  FORMAT! 'R=',F6.0, 'S') 

K S*0 

D05K=1,NHF,NHP 

KS=KS+1 

ENCODE (11  ,325  ,HLAB)HF  (K) 

5 CALL  LINES(HLAB,IPAK2,KS) 

325  F0RMAT('ALT»',F6.0,'S') 

001  OK  = 1 ,NR 
AVEC(K)*0 
CA2 (K )=C 
D 01  01  -1  ,NHF 
10  PR0B(K,I)=0 
D015L=1 ,NT 

READ ( 5, 104) ( A VE  W (K) , K *1 ,NR ) 

READt  5,104) (CA2SD (K ) ,K  = 1 , NR) 

104  F ORM A T ( 1 2 F6  • 0 ) 

READ(5,1GT)((PROBA(KtI),I*1,NHF)tK»1,NR) 

103  FORMAT  (12F4.3) 

501 5K -1 , NR 

AVEC(K)*AVEC(K)*AVEw(K) 

CA2(K)*CA2(K) *CA2SD(K) 

DOl  51  =1 , N H F 

15  PROB  (K,  I ) =PROB(K  , D + PPOBA  ( K,I  ) 
o o 2 Ok  = 1 ,np 

AVEC (K)SAVEC (K) /NT 
C A 2 ( K )*CA2(K) /NT 
D 0201  =1  ,NHF 

20  PROB < K, I ) =PR0B(K , I ) /NT 
WRITE (6,2C0)NTN 

200  FORMATtlHl'JVE  PROB  TABLE  ',A6,'  'EPRAIN ') 

C WRITE  AVE  PROB  TAPLE 
NHF 1 * 1NT ( NHF / 10) 

L2»0 

IF(NMF.LT.10)G0  TO  26 
D 035  I =1  , NMF  1 
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55 ' 

56  < 

57 

58 

59 
6C 
61 
62 

63 

64 

65 

66 

67 

68 

60 

70 

71 

72 

73 

74 

75 
76' 

77 

78 
70 
8 C 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
1C1 
1C2 
103 
1C4 

105 

106 
1C7 
1C8 
100 
110 
111 
112 
113 


Ll*L2*1 
L 2* 10  *1 

WRITE (6,206) (HF(L),L*L1 ,L2> 

206  F0R*AT(1MC9HALTITU0E*10F10.C) 

WRITE(6,237) 

207  F0RNAT(1hCShRANGE  ) 

D 03  5K  *1 , n r 

35  WRITE  (6,2  08)  R (K>  , (PROB  (K,J  >,J*Ll,L2> 

208  FORNAT(1HOF6.0,4X10F10.3> 

36  NHF2*NHF-1C»NHF1 
IF(NHF2.E9.0)GO  TO  39 
Ll*L2»1 

L 2*L2  +NH  F 2 

WRITE (6,206) (HF (L) , L = L 1 ,L2  ) 

WRITE  (6,207) 

0 03 6K  =1 ,NR 

38  WRITE  (6,2C8)R  (a)  , (PR  OP (X , J ),JSL1,L2) 

C PLOT  AVE  PROB  AS  FN  OF  ALTITUDE 

39  CALL  PH TS OR (2. ,5. 5) 

ENCODE(35 ,315 ,LA6ELP)NTN 

315  FORHAT('AVE  PROB  LOS  ',A6,'  TERRAIN','*') 

CALL  RESET('INTAXS') 

• CALL  XINTAX 

CALL  TITLE!'  ',1, 'ALTITUDE  ,HJ', 100, 'PROBABILITY*', ICQ, 3. 75, 2.75) 
CALL  NESSA6(LABELP,10C,-0.25,3.5) 

CALL  GRAF  (3. , 1000. , 5 0C0.  ,0  .,0.2  ,1  .0) 

0 C45K  =1  ,NR,NRP 
D044J  =1  , N HF 

44  p (j  >*pRce  (x , j ) 

45  CALL  CURVE(HF  ,P,NHF ,NHP) 

NRL«INT (NR/NRP) 

CALL  LEGE  NO  ( I PA  X 1 ,NRL, 4.0, 0.25) 

CALL  ENDGR(O) 

C PLOT  AVE  PROB  AS  FN  OF  RANGE 
CALL  PH TS OR (2., 1.5) 

CALL  TI TL  E ( ' ',1, 'RANGE,**', 100, 'PROBABILITY*', 100, 3.75, 2.75) 
CALL  GRAF (0. ,2000  .,16  300. ,0. ,0.2, 1.0  ) 

CALL  XINTAX 
D 055J  =1  , NHF  ,NHP 
D 05  OX  =1  ,NR 
50  P(K)»PR0B (X,J) 

55  CALL  CURVE (R ,P,NR ,NRP) 

NHL«I NT  ( N H F / N HP  ) 

CALL  LEGEN0(IPAK2, NHL, 4. 0,0.25) 

CALL  ENDGR(O) 

CALL  ENDPL(O) 

C PLOT  AVE  CRITICAL  ALT  AS  FN  OF  RANGE 
CALL  PH TS OR (2. ,5. 5) 

CALL  INTAXS 

CALL  TITLE!'  ', 1 , 'RANGE, NS ', 100 ,'CRI TICA L A LT I TU D E , N * ' , 1 CO , 3 . ’5 , 
X 7 5) 

C ALL  GRAF  (0  . , 200  0.,  16  300.  , C. , 2000  .,6 GOG.  ) 

CALL  CURVE(R ,AVEC  ,NR  ,NRP) 

CALL  CURVE(R ,CA2 .NR  .NRP) 

CALL  ENDGR(O) 

CALL  ENDPL(O) 

95  CONTINUE 

CALL  DONE  PL 
END 


■ 
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